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ABSTRACT 

One of the possible methods to distinguish among various dark matter candidates is 
to study the effects of dark matter decays. We consider four different dark matter 
candidates (light dark matter, gravitinos, neutralinos and sterile neutrinos), for each 
of them deriving the decaying/annihilation rate, the influence on reionization, matter 
temperature and CMB spectra. We find that light dark matter particles (1-10 MeV) 
and sterile neutrinos (2-8 keV) can be sources of partial early reionization (z ^ 100). 
However, their integrated contribution to Thomson optical depth is small ( ;$ 0.01) 
with respect to the three year WMAP results (r e = 0.09 ± 0.03). Finally, they can 
significantly affect the behavior of matter temperature. On the contrary, effects of 
heavy dark matter candidates (gravitinos and neutralinos) on reionization and heating 
are minimal. All the considered dark matter particles have completely negligible effects 
on the CMB spectra. 
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1 INTRODUCTION 

The nature of dark matter (DM) is one of the crucial open 
questions in cosmology. In the so called cold dark matter 
(CDM) theory DM particles are defined as 'cold' particles, 
because of their negligible free-streaming length (i.e. the 
length below which DM fluctuations are suppressed). The 
most famous alternative model to CDM is called warm dark 
matter (WDM), where DM particles are defined as 'warm' 
because of their longer free-streaming length. In WDM sce- 
narios the velocity dispersion of the particles is sufficient 
to smear out the fluctuations up to galactic scales, depend- 
ing on the mass of the particles (Padmanabhan 1995). This 
means that WDM models can alleviate the so called sub- 
structure crisis, which represents one of the most serious 
problems of CDM theories (Bode, Ostriker & Turok 2001; 
Ostriker & Steinhardt 2003). At present, there is no defini- 
tive evidence which allows us to exclude one of the two sce- 
narios and even the properties (mass, lifetime, etc) of cold 
and warm dark matter particles are substantially unknown. 

From an observational point of view, one of the most 
direct ways to detect DM particles and, maybe, distinguish 
between CDM and WDM is represented by particle decays 
(Chen & Kamionkowski 2004; Pierpaoli 2004). A small frac- 
tion of DM particles is expected to decay, the lifetime of 
this process generally depending on their density and mass. 
Many of the possible decay channels (Dolgov 2002) involve 



the emission of photons at wavelengths depending on the 
particle mass. So, in principle, it is possible to distinguish 
among various DM models, depending on the characteristics 
of emitted photons. It has also been pointed out (Sciama 
1982; Hansen & Haiman 2004; Chen & Kamionkowski 2004; 
Kasuya, Kawasaki & Sugiyama 2004; Kasuya & Kawasaki 
2004; Padmanabhan & Finkbeiner 2005; Mapelli & Ferrara 
2005; Zhang et al. 2006) that photons due to particle de- 
cays or annihilations can be sources of partial early reion- 
ization. In addition, the heating induced by particle decays 
can leave an imprint in the 21 cm background, which could 
be detected by new generation of radio telescopes (LOFAR, 
PAST, SKA, etc.). 

At the moment, constraints on the radiation emitted 
by particle decays are loose. The SPI spectrometer aboard 
ESA's INTEGRAL satellite recently detected an excess in 
the 511 keV line emission, due to positron-electron annihi- 
lation, from the galactic bulge (Knodlseder et al. 2005). The 
only two viable explanations are that this excess is due to 
positrons produced by thermonuclear Type la supernovae 
(Dermer & Murphy 2001) or by decaying/annihilating dark 
matter (Hooper & Wang 2004; Boehm et al. 2004). This last 
hypothesis, although exotic, has triggered many theoretical 
studies (Ascasibar et al. 2006; Kawasaki & Yanagida 2005; 
Kasuya & Takahashi 2005; Casse & Fayet 2005; Kasuya & 
Kawasaki 2006), which are aimed to put constraints on var- 
ious DM properties by using SPI/INTEGRAL observations. 
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In this paper, we consider some of the most popular 
cold and warm DM particles, calculating their approximate 
decaying rate (Section 2), their influence on the ionization 
fraction, on the Thomson optical depth, on the behavior of 
matter temperature (Sections 3 and 4) and on the cosmic 
microwave background (CMB) spectra (Section 5). 
In all the cases, we make the assumption that the DM is 
composed of one single species of particles. We consider only 
'standard' DM candidates, neglecting more exotic scenarios 
(such as Q-balls, light scalar bosons, etc). In particular, we 
study three different candidates for the case of CDM: (i) the 
axino, as representative of light dark matter (LDM; 1-100 
MeV; Hooper & Wang 2004), (ii) the gravitino and the (iii) 
neutralino, as heavy dark matter candidates ( £ 100 MeV). 
For the WDM we consider only the sterile neutrino, reduc- 
ing our analysis to its radiative decay channel. In fact, we 
do not pretend to present a complete overview of DM can- 
didates. Instead, we would like to give a basic description of 
the effects of DM decays, taking as an example some of the 
standard DM candidates. Our aim is to point out the differ- 
ences among the considered DM particles, with particular 
care for cosmic reionization and heating. 



2 METHOD 

For each considered particle model we derived the energy 
injection rate per hydrogen nucleus, (dm, as follows. 

In comoving coordinates the photon emission rate due 
to particle decay can be generally written as: 

^ _ H2. J*(o)-t(z)]/T 
At T 



(1) 



where no and r are the current density and the lifetime of 
particles, respectively, and t{z) ~ | H^ 1 ^^ 2 {1 + z)~ 3 / 2 
(neglecting Qa) is the time elapsed from the Big Bang to 
redshift z. 

Instead, in the case of annihilations, the photon emis- 
sion rate in comoving coordinates is 



— = n (1 + z) avC, 
At 



(2) 



where a v is the annihilation cross-section (see Section 3.3 
for details) and C is the clumping factor. To quantify C is 
beyond the scope of this paper. However, especially at very 
high redshift (z 3> 10), where annihilations play their most 
important role, we can roughly approximate C ~ 1. 
Then, in both the cases (dm is simply: 



An E. 



iDM = 



At n? (3) 

where _E 7 ~ uidm/2 is the energy of the emitted photon 
(rriDM being the mass energy of the DM particle), and nt 
the current density of baryons (we take n& = 2.7 x 10~ 7 
cm~ 3 , Spergel et al. 2003). 

The energy (dm partially goes into ionizations of hydro- 
gen and helium atoms and partially into heating. We adopt 
the rough approximation by Chen & Kamionkowski (2004), 
for which a fraction (1 — a;)/3 of €dm contributes to ioniza- 
tions and a fraction (l + 2cc)/3 goes into heating {x being the 
ionized fraction). This approximation is correct if the DM 
decays/annihilations produce (directly or via secondary in- 
teractions) photons whose energy is mostly absorbed by the 



intergalactic medium in a Hubble time. As fig. 2 of Chen 
& Kamionkowski (2004) shows, photons emitted by sterile 
neutrinos (~ 1 — 4 keV) fully satisfy this requirement. Parti- 
cle decays/annihilations which produce electrons of energy 
lower than ~ 1 GeV (Chen & Kamionkowski 2004) equally 
satisfy this requirement, because these electrons can lead (by 
collisional ionization, excitation or inverse-Compton scatter- 
ing) to photons which are immediately absorbed. Then, also 
for LDM and gravitinos we can adopt the approximation of 
Chen & Kamionkowski (2004). For neutralinos (whose mass 
is higher than 30 GeV) this approximation is far too opti- 
mistic; however we can consider it as an upper limit. In prac- 
tice, in the optically thin case most of the decay /annihilation 
products can propagate to redshift and can be constrained 
by observing the X-ray /gamma-ray background (Chen & 
Kamionkowski 2004). 

We made our calculations using an upgraded version of 
the public code RECFAST (Seager, Sasselov & Scott 1999, 
2000). In particular, we modified the evolution equations as 
follows (Padmanabhan & Finkbeiner 2005). 
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3(1 + /He) 

where xh (xHe) is the ionized fraction of hydrogen (helium) 
atoms, Eth, h = 13.6 eV {Eth, He = 24.6 eV) is the ionization 
energy of hydrogen (helium) atoms, /n e is the helium-to- 
hydrogen ratio by number, Tm is the matter temperature, 
k,B the Boltzmann constant and £ = [H{z){l + z)] _1 - 



3 COLD DARK MATTER 

First we consider CDM particles. Heavy CDM particles (> 
100 MeV) are not considered a viable source for the 511 keV 
emission in the galactic center. In fact, in the case of neu- 
tralinos (with mass higher than 30 GeV), the request of a 
sizable R-parity violation, needed to allow considerable neu- 
tralino decays, would determine a too short lifetime and the 
neutralino would cease to be a good DM candidate (Hooper 
& Wang 2004). On the other hand, gravitino decays are pos- 
sible; but the gravitino lifetime is far too long to match the 
511 keV emission from the galactic center (Hooper & Wang 
2004). Then, previous studies (Bohem et al. 2004; Hooper 
& Wang 2004; Ascasibar et al. 2006) proposed light cold 
DM candidates (1-100 MeV) to be sources of the 511 keV 
emission from the galactic center. In the following, we will 
consider firstly LDM particles (axinos) and secondly heavy 
DM particles (gravitinos and neutralinos). 



3.1 Light dark matter 

LDM candidates (1-100 MeV) can produce positrons both 
via decay (axinos, Hooper & Wang 2004) and via annihila- 
tion (Bohem et al. 2004). They can easily explain the 511 
keV emission from the galactic center and satisfy the DM 
relic density {Qdm ~ 0.23, Spergel et al. 2003). The upper 
limit of their mass, if they are the source of the 511 keV line, 
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is probably much less than 100 MeV, due to constraints on 
the bremsstrahlung emission (Beacom, Bell & Bertone 2005; 
Casse & Fayet 2005; Beacom & Yiiksel 2006). 

Hooper & Wang (2004) derived in a very simple way the 
lifetime of decaying LDM particles (i.e. axinos) necessary to 
produce the observed 511 keV emission: 



, , 1n 26 / ttildm 
4 x 10 s 



/ ttildm \ 
\ MeV J 



(7) 



where ttibm is the mass of a LDM particle. Under our as- 
sumptions, the current density of LDM particles can easily 
been derived as: 



n 



Pc 



TtlLDM 



(8) 



where p c is the critical density of the Universe. Substituting 
equation (7) and (8) into equation (1) and implementing it in 
RECFAST (through equations 4-6), we derive the influence 
of LDM particles on ionization and heating, shown in Fig. 1. 
The contribution of LDM starts to be important at redshift 
z < 100. The current value of x e should be of the order of 0.1 
(for masses ttildm ^ 5 MeV). Then LDM does not produce 
a complete reionization; but can be an important source of 
early partial ionization. Also the matter temperature Tm 
starts to differ from the case without LDM at z ~ 100. 
At z ~ 20 T M is ~ 200 K, a factor ~ 20 higher than in 
the unperturbed case. At lower redshift Tm becomes of the 
order of 10 4 K. These high temperatures enhance the role of 
collisions (both between electrons and hydrogen atoms and 
between two hydrogen atoms). However, the collision time- 
scale (Palla, Salpeter & Stahler 1983) remains always two 
or more orders of magnitude longer than the Hubble time, 
allowing us to neglect collisions in our calculations. 



3.2 Gravitinos 

The most probable gravitino masses are m 3/ / 2 < IkeV and 
m 3/2 > ITeV (Nowakowski & Rindani 1995). In the first 
case, the gravitinos are the lightest supersymmetric parti- 
cle (LSP) and then are stable: they are good WDM can- 
didates; but their decay rate is negligible. In the second 
one, gravitinos are so unstable that they decay in the early 
universe, and, as a consequence, they are not viable DM 
candidates. Some models assume that gravitinos have mass 
m 3/2 ~ 10 — 100 MeV, are the LSP and can violate the 
R-parity (Hooper & Wang 2004). In this case gravitinos are 
good DM candidates and have non-negligible decay rate. For 
this case, we calculated the gravitino contribution to reion- 
ization and heating, assuming lifetime: 



10 31 s 



100 GeV 



0.1 GeV 

m-3/2 



(9) 



where mr is the slepton mass and A is the R-parity violating 
leptonic trilinear coupling. The current density of graviti- 
nos can be derived as indicated for LDM. We find that, 
because of such a long lifetime, the contribution of graviti- 
nos to heating and reionization is negligible (Fig. 2) . For the 
same reason, Hooper & Wang (2004) show that gravitinos 
are unable to produce the 511 keV excess from the galactic 
center. 
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Figure 1. Ionized fraction (bottom panel), Thomson optical 
depth (central panel) and matter temperature (upper panel) as 
a function of redshift due to decaying LDM of masses 1 (thick 
dotted line), 5 (dashed) and 10 MeV (solid). The thin solid line 
represents, from bottom to top, the relic fraction of free electrons, 
their contribution to Thomson optical depth and the matter tem- 
perature without particle decays. 
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Figure 2. Ionized fraction (bottom panel), Thomson optical 
depth (central panel) and matter temperature (upper panel) as 
a function of redshift due to decaying gravitinos of masses 10 
(dashed) and 100 MeV (solid). The thin solid line is the same as 
in Fie. 1. 
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3.3 Heavy dark matter: neutralinos 



Here we will consider as 'heavy' dark matter the neutralinos. 
It is a merely indicative classification, given the uncertainties 
on the various models. The discussion of the details of dif- 
ferent supersymmetric models is beyond the purpose of this 
paper. Neutralinos are thought to be very massive (m x > 30 
GeV). So, if they could decay (violating the R-parity), their 
lifetime should be very short, and they could not be a vi- 
able DM candidate. Then, the neutralino, if exists, must 
be perfectly stable, and we will not treat neutralino decay. 
However, neutralinos can annihilate. The annihilation cross- 
section is generally fit by (Bertone, Hooper & Silk 2005): 

av = a + bv 2 + 0{v 4 ), (10) 

where a and b are constant, whose values are constrained 
by the DM relic density condition, and v is the neutralino 
velocity, which depends on the DM temperature and thus 
on the redshift. In the present epoch neutralinos are non- 
relativistic, then the current annihilation cross-section can 
be written as a v ~ a. However, the cross-section at the 
freeze-out time should depend on v and be higher than the 
current value. As a rough approximation, Padmanabhan & 
Finkbeiner (2005) consider a thermally averaged, redshift 
independent cross-section, (uv) = 2 x 10 -26 cm 3 s _1 . For 
comparison, we made the same assumption. Then, the an- 
nihilation rate becomes: 

? = 2-< 
dt 



x 10- 42 (1 



n3 -3 

z ) cm s 



n 



2 x 10" 8 c 
(av) 



2 x 10" 26 cm 3 s- 



(11) 



where no 



p c /m x . In our calculations we assume 



m x =100 GeV. We have implemented this equation into 
RECFAST. The contribution of neutralino annihilations 
both to ionizations and heating is negligible (Fig. 3; dashed 
line) for (av) — 2 x 10 -26 cm 3 s" 1 (in agreement with Pad- 
manabhan & Finkbeiner 2005). As an upper limit, we con- 
sidered also the case where (uv) = 10 -24 cm 3 s _1 , which is 
the highest value to be consistent with the first year WMAP 
data (see Colafrancesco, Profumo & Ullio 2005) . Also in this 
case, the contribution to heating is negligible and the ion- 
ization fraction remains of the order of 10 -3 . However, an- 
nihilations are particularly important at very high redshift 
(z ^ 100), where the particle density is very high. For this 
reason, even if the ionization fraction due to annihilations 
remains always very low, the Thomson optical depth is sig- 
nificantly high (r e ~ 0.05), even more than for LDM. 

These results must be considered very optimistic upper 
limits. In fact we are assuming that nearly all the energy 
of the DM particle is immediately deposited into ionization 
or heating; whereas we expect that the electrons produced 
by neutralino annihilations Compton-scatter the CMB pho- 
tons up to a energy ~ 1 — 10 (1 + z) MeV, which cannot be 
significantly absorbed by the intergalactic medium within a 
Hubble time (Chen & Kamionkowski 2004). 



WARM DARK MATTER: 
NEUTRINOS 



STERILE 



Sterile neutrinos are one of the most popular WDM candi- 
dates (Colombi, Dodelson & Widrow 1996; Sommer-Larsen 
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Figure 3. Ionized fraction (bottom panel), Thomson optical 
depth (central panel) and matter temperature (upper panel) as 
a function of redshift due to neutralinos for (a v) = 2 X 10 — 26 
(thick dashed line) and 10 — 24 cm 3 s — 1 (solid). In both the cases 
the neutralino mass is 100 GeV. The thin solid line is the same 
as in Fie - . 1. 



& Dolgov 2001), even if Seljak et al. (2006) seem to exclude 
that they are the only component of DM on the basis of Ly- 
a forest power spectrum measurements. They can exist only 
if neutrinos have non-zero mass and mixing angles, as pre- 
dicted by the standard oscillation theory (Dolgov & Hansen 
2002; see Dolgov 2002 for a complete review of sterile neu- 
trino properties). There are many possible decay channels 
of sterile neutrinos (Dolgov 2002). In this paper we are in- 
terested on the radiative decay, i.e. the decay of a sterile 
neutrino into a lighter neutral fermion (such as an active 
neutrino) and a photon, because of its effects on the cos- 
mic ionization and heating. From the comparison between 
the predicted background flux due to radiatively decaying 
sterile neutrinos and the hard X-ray background (Bauer et 
al. 2005), Mapelli & Ferrara (2005) have established an up- 
per limit of 14 keV for the sterile neutrino mass. This limit 
can now be lowered to m v „ < 11 keV, adopting the relation 
between the mixing angle and the mass recently derived by 
Abazajian (2006). A stronger upper limit, m v , < 8.2 keV, 
has been derived from X-ray observations of the Virgo clus- 
ter (Abazajian 2006). Furthermore, Viel et al. (2005) derived 
a lower limit m„ s > 2 keV from the study of matter power 
spectrum fluctuations. Then, sterile neutrino masses are al- 
lowed from 2 to 8 keV, a very narrow range. 

The lifetime for sterile neutrino radiative decay is 
(Mapelli & Ferrara 2005): 



512 7T 4 
9 ^em 



G F ml 



(12) 



where a em is the fine structure constant, Gf the Fermi con- 
stant, ra„ s the sterile neutrino mass and sin# the mixing 
angle. To derive sin 6 we adopt the following relation (Abaza- 
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Figure 4. Ionized fraction (bottom panel), Thomson optical 
depth (central panel) and matter temperature (upper panel) as a 
function of redshift due to radiatively decaying sterile neutrinos 
of masses 2 (thick dotted line), 4 (dashed) and 8 keV (solid). The 
thin solid line is the same as in Fig. 1. 



jian 2006): 



x <^ 0.527 erfc 



/ 3-4keV \ / Qdm \ 
\ m va ) \ 0.26 J 

( Tqqd \ 21i 
V 170 MeV 7 



1/2 



-1.15 



(13) 



where Qdm is the dark matter density and Tqcd the tem- 
perature of quark-hadron transition. 

Assuming that all the DM is composed by sterile neu- 
trinos and substituting equation (12) into equation (1) we 
derive through RECFAST the ionization and heating his- 
tory also for WDM particles (Fig. 4). Also sterile neutrinos 
start to play a role into the reionization and heating at red- 
shift 2 ~ 100, and their behavior is close (even if the global 
contribution is slightly lower) to that of LDM. 



5 EFFECTS ON THE CMB SPECTRUM 

In the previous section we have shown that decaying DM, 
and especially LDM and sterile neutrinos, can modify the 
ionization fraction, with respect to the value due to relic 
electrons, already at high redshift. This fact should leave 
some imprint on the CMB spectrum (Chen & Kamionkowski 
2004; Pierpaoli 2004; Padmanabhan & Finkbeiner 2005). To 
check whether these effects are measurable, we simulated the 
expected CMB spectrum in the case we take into account 
DM decays. This has been done by implementing our mod- 
ified version of RECFAST in the version 4.5.1 of the public 
code CMBFAST (Seljak & Zaldarriaga 1996; Seljak et al. 
2003). 

Fig. 5 shows the temperature - temperature (TT), tem- 
perature - polarization (TE) and polarization - polariza- 
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Figure 5. Temperature-temperature (top panel), polarization- 
polarization (central panel) and temperature-polarization (bot- 
tom panel) spectra. Thick lines indicate the CMB spectrum de- 
rived assuming Thomson optical depth r e = 0.09 and a sudden 
reionization model (consistent with the three year WMAP data); 
thin lines indicate the CMB spectrum derived assuming r e = 0. 
Dashed (solid) lines indicate the CMB spectrum obtained (with- 
out) taking into account the decays of 10 MeV LDM particles. 
The two thick lines, solid and dashed, appear superimposed, be- 
cause the contribution of decaying particles (the dashed line) is 
completely hidden by the stronger effect of a sudden reionization 
with r e = 0.09. Open circles in all the panels indicate the three 
year WMAP data (Hinshaw et al. 2006; Page et al. 2006; Spergel 
et al. 2006). 



tion (EE) spectra, in the case of 10 MeV decaying LDM (i.e. 
the particle for which we achieved the maximum contribu- 
tion to the reionization among the considered ones), com- 
pared with the recent three year WMAP data (Spergel et al. 
2006; Page et al. 2006). The contribution due to DM decays 
alone is negligible. There is a sensible difference only in the 
lowest multipoles (I < 10) of the EE spectrum. This effect 
can be seen in the central panel of Fig. 5, where the thin lines 
show the expected EE spectra by considering (dashed line) 
and neglecting (solid line) DM decays, respectively. This ef- 
fect, small and concentrated at low multipoles, is justified 
by the fact that DM decays produce a very small Thom- 
son optical depth (r e ^ 0.01) and that they are important 
especially at very low redshift, due to their long lifetime. Is 
such a modification of the EE spectrum measurable? If there 
are other sources of reionization besides DM decays (as it 
seems to be likely, considering the Thomson optical depth, 
r e = 0.09±£;^, measured by WMAP; Spergel et al. 2006), 
the influence exerted on the EE spectrum by the decaying 
DM would be completely hidden by the stronger effects due 
to these other reionizing sources. This can be seen in Fig. 5, 
where the thick lines show the TT/EE and TE spectra as- 
suming r e = 0.09 in the case with (dashed line) and without 
(solid line) DM decays. We found that the effects of DM de- 
cays are washed out by those of other reionizing sources also 
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for lower values of r e consistent with the three year WMAP 
results (down to r e = 0.06, corresponding to a sudden reion- 
ization at z ~ 6). 

Because 10 MeV LDM particles produce the highest 
ionization fraction among the considered models, the effects 
on the CMB spectra due to other species of DM particles 
will be far more negligible. However, stronger effects on the 
CMB spectra can be due to annihilations of LDM particles 
(Zhang et al. 2006), which are not considered by this paper. 



6 CONCLUSIONS 

We examined the contribution to cosmic reionization and 
heating of different models of decaying/annihilating DM. In 
the case of quite light particles ( 10 MeV) this contribu- 
tion is significant. Light particles (LDM, or sterile neutri- 
nos) become important at z ~ 100, as they provide an early 
source of (partial) ionization for the intergalactic medium. 
They are expected to produce a Thomson optical depth 
r e < 0.01 (r e < 0.001) in the case of LDM particles (ster- 
ile neutrinos), which is smaller than the value derived from 
the three year WMAP data (r e = 0.09), but non- negligible. 
Changes in the matter temperature are also important, and 
the role of DM decays on the history of 21 cm emission 
should be investigated. On the contrary, heavier particles 
(gravitinos and neutralinos) do not have significant influ- 
ence on reionization and heating. 

This result could be crucial in distinguishing between 
light and heavy dark matter models, if new measures will be 
available of the reionization history and/or the behavior of 
the matter temperature (e.g. mapping the 21 cm emission at 
z ~ 10 — 50). However, it is quite impossible to distinguish 
among different species of light DM particles, such as sterile 
neutrinos or LDM. 

Finally, in the case of light particles, an early (z ^> 
20) increasing of the ionization fraction and of the baryon 
temperature could catalyse the production of H2 and HD 
molecules, affecting the entire history of structure formation 
(Shchekinov & Vasiliev 2004; Biermann & Kusenko 2006). 
On the contrary, no constraints on DM particles can be de- 
rived from their effects on the CMB spectra. 
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